quad('x.*log(1+x)',0,1)